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Characterization of majorization monotone quantum 

dynamics 

Haidong Yuan 


Abstract —In this article I study the dynamics of open quantum 
system in Markovian environment. I give necessary and sufficient 
conditions for such dynamics to he majorization monotone, which 
are those dynamics always mixing the states. 


I. Introduction 

In the last two decades, control theory has been applied 
to an increasingly wide number of problems in physics 
and chemistry whose dynamics are governed by the time- 
dependent Schrodinger equation (TDSE), including control of 
chemical reactions Cl, m, 10, El, 0, 0, Q, 0, state- 
to-state population ttansfer 191, IfTOl , IfTTI . IfT^ . ||T31 , shaped 
wavepackets lfT4l . NMR spin dynamics ITSl , ifTfil , iflTl , flSl , 
CSl, Bose-Einstein condensation iMl, El, El, quantum 
computing 1231 . l24l . l25l . l2^ . l27l . oriented rotational 
wavepackets l28l . l29l . l30l . etc. More recently, there has 
been vigorous effort in studying the control of open quantum 
systems which are governed by Lindblad equations, where 
the central object is the density matrix, rather than the wave 
function El, El, IESI, El, Ea, l36l, El. The Lindblad 
equation is an extension of the TDSE that allows for the 
inclusion of dissipative processes. In this article, I will study 
those dynamics governed by Lindblad equations and give 
necessary and sufficient conditions for the dynamics to be 
majorization monotone, which are those dynamics always 
mixing the states. This study suggests that majorization may 
serve as time arrow under these dynamics in analog to entropy 
in second law of thermal dynamics. 

The article is organized as following: section|II]gives a brief 
introduction to majorization; section |III] gives the definition of 
majorization monotone quantum dynamics; then in section HVl 
necessary and sufficient conditions for majorization monotone 
quantum dynamics are given. 

II. Brief introduction to majorization 

In this section I give a brief introduction on majorization, 
most stuff in this section can be found in the second chapter 
of Bhatia’s book l42l . 

Eor a vector x = (xi, ...,Xn)'^ in K.", we denote by x^ = 
(a:f, a permutation of x so that xf > xj if i < j, 

where 1 < i, j < n. 
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Definition 1 (majorization) A vector x G R” is majorized 
by a vector y G K” (denoted by x ^ y), if 

( 1 ) 

i=i 

for d = 1,... ,n — 1, and the inequality holds with equality 
when d = n. 

Proposition 1 x -< y iff x lies in the convex hull of all Piy, 
where Pi are permutation matrices. 

Proposition 2 x ^ y if and only if x = Dy where D is 
doubly stochastic matrix. 

Remark 1 A doubly stochastic matrix D is a matrix with non¬ 
negative entries and every column and row sum to 1, i.e., 

dij > 0 , dij = 1 , dij = 1 . 

Proposition 3 Suppose / is a convex function on R, and x -< 
y in R”, then 

n n 

i i 

Proposition 4 Eor a vector A = (Ai,..., A„)^, denote D\ 
a diagonal matrix with (Ai,...,A„) as its diagonal entries, 
let a = (ai,...,a„)^ be the diagonal entries of matrix 
A = D\K, where K G SO(n). Then a -< A. Conversely 
for any vector a ^ A, there exists a K G SO{n), such that 
(oi, ...jOn)^ are the diagonal entries of A = K^D\K. 

Remark 2 SO(n) is the group of special orthogonal matrices, 
K G SO{n) means K^K = I and det{K) = 1. 

III. Majorization in open quantum dynamics 

The state of an open quantum system of N-level can be 
represented by a x positive semi-definite, trace 1 matrix, 
called density matrix. Let p denote the density matrix of an 
quantum system, its dynamics in markovian environment is 
governed by the Lindblad equation, which takes the form 

p=-i[H{t),p]+L{p), (2) 

where —i[H, p] is the unitary evolution of the quantum system 
and L{p) is the dissipative part of the evolution. The term L{p) 
is linear in p and is given by the Lindblad form ESI, ESI, 

L{p) = '^ac,0{FcpFl - \{FlFa,p}), 

0.0 
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where Fa, Fp are the Lindblad operators, which form a basis 
of N X N trace 0 matrices (we have N'^ — 1 of them) and 
{A, B} = AB + BA. If we put the coefficient a^/S into a 
(iV^ — 1) X (iV^ — 1) matrix A = ( 0 ^/ 3 ), it is known as the 
GKS (Gorini, Kossakowski and Sudarshan) matrix which 
needs to be positive semi-definite. 

Eq. ^ has the following three well known properties: 1) 
Tr(p) remains unity for all time, 2) p remains a Hermitian 
matrix, and 3) p stays positive semi-definite, i.e. that p never 
develops non-negative eigenvalues. 

Definition 2 Suppose pi and p2 are two states of a quantum 
system, we say pi is majorized by p 2 {pi -< P 2 ) if the 
eigenvalues of pi is majorized by the eigenvalues of p2 
(A(pi) ^ \{P 2 ))- 

Basically majorization gives an order of mixed-ness of 
quantum states, i.e., if pi -< p 2 , then pi is more mixed than 
P2, which can be seen from the following propositions. 

Definition 3 (Von Neumann entropy) The Von Neumann 
entropy of a density matrix is given by 

S{p) = -Tr[plog{p)]. 

Proposition 5 If pi -< p 2 , then S{pi) > S{p 2 ). 

Proposition 6 If pi F p 2 , then Tr{p\) <Tr{p 2 ). 

Remark 3 The above two propositions can be easily derived 
from Proposition [3] 

The entropy and trace norm are usually used to quantify how 
mixed quantum states are. But majorization is a more strong 
condition than these two functions, and in some sense it gives 
a more proper order of mixed-ness as we can see from the 
following proposition. 

Proposition 7 l\44\l pi -< p 2 if and only if pi can be 

obtained by mixing the unitary conjugations of p 2 , i.e., pi = 
YjiPiUiP 2 U}, where pi > 0, = 1 and Ui are unitary 

operators. 

IV. Necessary and sufeicient condition of 

MAJORIZATION MONOTONE QUANTUM DYNAMICS 

Definition 4 (Majorization monotone dynamics) An open 
quantum dynamics governed by Eq. © is majorization mono¬ 
tone if and only if p{t 2 ) F p{ti) when t 2 > ti, 

Intuitively majorization monotone dynamics are those kind of 
dynamics which always mixing the states. As we can see from 
Proposition |5] these kind of dynamics always increase the 
entropy of the system. One can immediately see a necessary 
condition for a dynamics to be majorization monotone: the 
state Pi — has to be a steady state of such dynamics, 
where I is identity matrix. As AI is the most mixed state. 


any state pi < Aj would imply pi = ^7. The question now 
is whether this condition is also sufficient. 

Let’s first look at a simple system: a single spin in marko- 
vian environment. 

A. An example on single spin 

Take the general expression of the master equation 

p = -i[H,p\+ L{p), ( 3 ) 

where 

^(p) = '^aap{FapFl - \{FlFa,p}). 

OL0 

Eor the single spin, we can take the basis {Fa} as normalized 
Pauli spin operators iTy,Uz}, where ax = (i q), ay = 

(° “*), and CTz = (J _?]^). The coefficient matrix 

( ^xx ^xy 
Qiyx ^yy 
^zx ^zy 

is positive semi-definite. 

If identity state is a steady state, the right hand side of Eq. 
© should be 0 when p = ^I. As —i[H, ^/] = 0, so the 
condition reduces to L{I) = 0 (since L{p) is a linear map, 
we can ignore the constant A). This is equivalent to 

Y,aap[Fa,Fl]=f). ( 4 ) 

a,/3 

In the single spin case, substitute F by Pauli matrices and it 
is easy to see that the above condition reduces to 

i.e., the GKS matrix should be real symmetric, positive semi- 
definite matrix BTll . while the general GKS matrix is Her¬ 
mitian, positive semi-definite. We want to see whether the 
dynamics of single spin under this condition is majorization 
monotone. 

As majorization monotone is defined by the eigenvalues of 
density matrices, we are going to focus on the dynamics of the 
eigenvalues of density matrix. Let A be its associated diagonal 
form of eigenvalues of density matrix p, i.e., A is a diagonal 
matrix with eigenvalues of p as its diagonal entries. At each 
instant of time, we can diagonalize the density matrix p{t) = 
U{t)A{t)U^(t) by a unitary matrix (7(f). 

Substitute p(t) = U{t)A{t)U'^(t) into Eq. ©, we get 

pit) =U{t)A{t)U^{t) + Uit)A{t)U\t) 

+ Uit)Ait)U^{t) 

= - iH'it)Uit)A{t)U\t) + Uit)A{t)U\t) 

+ Uit)Ait)U\t)iH'{t) 

= - i[H'it), U{t)A{t)U\t)\ + U{t)A{t)U^{t) 

= - i[Hit),Uit)A{t)U\t)\ + L[(7(f)A(f)(7t(f)], 
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where H'{t) is defined by (j{t) = —iH'(t)U(t), which is 
Hermitian. We obtain 

A(t) =uHt){-i[H{t) - H'it),U{t)A{t)U\t)] 

+L[umm\muit) 

+ U\t)L[U{t)A{t)U\t)]U{t). 

Note that the left side of the above equation is a diagonal 
matrix, so for the right side we only need to keep the diagonal 
part. It is easy to see that the diagonal part is zero for the first 
term, thus we get 

A{t) = diag{U'>{t)L[U{t)A{t)U'<{t)]U{t)), (7) 

where we use diag{M) to denote a diagonal matrix whose 
diagonal entries are the same as matrix M. 


A(t) = diag{U'<L{UAU'<)U) 

= diag(^^aap(U^ FqlJ AU^ FgU 
a/3 

-\{U^FIUU^F^U,A})) (8) 

= diagC^ao,p{U'< FaU AU"^ FpU 

a/3 

- i{C/tF^C/C/tF„C/,A})). 

For the last step we just used the fact that Fp is a Pauli matrix 
which is Hermitian. Now 


U'^F^U = c^^F^, 


where C = 


^xz 

Cyz 

Czz 


G S'0(3) is the adjoint 


'-'yx ^yy 
^zx ^zy 

representation of U. Substituting these expressions into Eq. 
we obtain 


A{t) = diagC^a'^p{Fc,AFp - -{FpFo,,A})), 


(9) 


a/3 

where 

are entries of transformed GKS matrix. 


A' = C^AC. 

We can write A{t) = ^ ^ = 1/ + X{t)az, where 

A(f) G [0, A]. Substitute it into Eq. ®, we obtain the dynamics 
for A(t), 

\{t) = -(a'n +a22)A(f). 

Erom Proposition |4] we know 


< «!! + 022 ^ fJ-2 + g-1-, 

where gi > ^2 ^ /^3 are eigenvalues of the GKS matrix. 
Prom this it is easy to see that at time T, the value of A(T) 
lies in the following interval 

[e-(^i+^=)^A(0),e"('^^+'"3)'^A(0)], 


which is always less or equal to A(0). And this is sufficient 
for the dynamics to be majorization monotone in the single 
spin case, as (A + A(r), A — X{T)) -< (^ + A(0), A - A(0)) 
when A(T) < A(0). So L{I) = 0 is necessary and sufficient 
for the dynamics to be majorization monotone in the single 
spin case. And this condition holds even we have coherent 
control on the spin, as the identity state remains as steady state 
in the presence of control as we can see from the following 
controlled dynamics: 

P = + ^u,H„p] + L{p), 

i 

where UiFli are our coherent controls. Suppose the con¬ 
trollers are able to generate any unitary operations on the spin 
fast compare to the dissipative rate, then A(T) can actually take 
any value in the interval 

i.e., the reachable set for the single spin under the controlled 
Lindblad dynamics is 

= (j.) ) [;t|A(T) e 

[e-('^i+'"^)'^A(0), U G SU{2)}. 

This is to say that although controls can’t reverse the 
direction of mixing, it can change the rate within some region. 

Remark 4 Prom the single spin case, we can see that a dy¬ 
namical system being majorization monotone does not imply 
the states of the system always converge to identity state, as 
identity being a steady state does not exclude the possible 
existence of other steady states, for example, the dynamics 
given by 

P= -i{crz,p] +-i[az,[az,p]], 

which describes the transverse relaxation mechanism in NMR, 
satisfies the majorization monotone condition in the single spin 
case, and it is easy to see that the state of this system does not 
necessary converge to identity matrix, in fact it can well be 
converged to any state of the form Aj-l-acr^, where a G [0, A]. 


B. General Case 

In this section, we will show that L{I) = 0 is also sufficient 
for the dynamics to be majorization monotone in the general 
case. Suppose we solved the Lindblad equation 

p=-i[H{t),p]+L{p), (10) 

integrated this equation from ti to t 2 , where t 2 > ti, and get 
a map: 

^ : p{ti) -)> p{t2). 

Such a map has a Kraus operator sum representation ll40l . 

HI: 

p{t2) = 'S{p{tl))=J2K^p{tl)Kl (11) 

i 

where in our case {Ki} are N x N matrices, which depend 
on the dynamical Eq. (fTOl i and the time difference between ti 
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and ^ 2 - Also the Kraus operator sum has to be trace preserving 
as the trace of density matrix is always 1, which implies that 

Y,kIk, = i. 

i 

If we have additional condition that identity state is a steady 
state of this dynamics, which means if p{ti) = then p{t 2 ) 
remains at substitute them into the Kraus operator sum 
representation, we will get an extra condition 

Y,k,kI = i. 

i 

We will show these two conditions are enough to ensure the 
dynamics to be majorization monotone. First let’s diagonalize 
p{ti) and p(f 2 ): 

p{ti) = UiK{p[ti))Ul, 

P{t2) = U2A{p{t2))Ul 


where A(p) are diagonal matrix with eigenvalues of p as its 
diagonal entries, substitute them into Eq. EJ, we get 

U2A{p{t2))Ul =Y,K,UMp{h))ulKl 

* ( 12 ) 

A(p(i2)) = ^ UlK,U^K{p{h))UlKlU2. 

i 

Let Vi = U^KiUi, then 

A{p{t2)) = Y.^^Hpitl))V^, 

i 

and it is easy to check that 

= uI{Y,k,k})U2 = I, 

i i 

Y.K^V = UliY,K,Kj)U,=L 


(13) 


(14) 


It is a linear map from the eigenvalues of p{ti) to eigen¬ 
values of p{t 2 ), so we can find a matrix D, such that 

\{p{t2))=D\{p{h)), (15) 

where A(p) is a vector in with eigenvalues of p as its 
entries, which is arranged in the same order as the diagonal 
entries of A{p). The matrix D can be computed from Eq. ( fT3l i: 

Dap =''^^\{Vi)al3\ , 
i 

where Dap and (Ei)a /3 are the a/3 entry of D and Vi 
respectively. It is straightforward to show that, by using the 
two conditions in Eq. (fT4l l. 

Dap = 1, 

a. 

Dap = 1, 

/3 

i.e., D is a doubly stochastic matrix. Erom Proposition |2l we 
get 

\{p{t2)) -< A(p(fi)), 


so 

p{t2) A p{ti)yt2 > fl, 

i.e. it is majorization monotone. 

Erom Proposition |5] and |6] it is easy to see that majoriza¬ 
tion monotone implies entropy monotone and trace norm 
monotone, and they share the same necessary and sufficient 
condition; L{I) = 0. 


V. Conclusion 

Understanding open quantum systems is an important prob¬ 
lem for a wide variety of physics, chemistry, and engineering 
applications. This paper analyzed the dynamics of open quan¬ 
tum systems and gives necessary and sufficient condition on 
majorization monotone dynamics, which are those dynamics 
always mixing the states. This suggests that for this class of 
dynamics, majorization defines an evolution arrow, which begs 
for the connection to the entropy arrow in the second law of 
thermal dynamics. I hope further investigation will reveal more 
on this connection. 
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